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ABSTRACT 

We present results from the application of a global photometric calibration (GPC) 
procedm'e to calibration data from the first 2 years of The Two Micron All Sky 
Survey (2MASS). The GPC algorithm uses photometry of both primary standards 
and moderately bright 'tracer' stars in 35 2MASS calibration fields. During the first 
two years of the Survey, each standard was observed on approximately 50 nights, 
with about 900 individual measurements. Based on the photometry of primary 
standard stars and secondary tracer stars and under the assumption that the nightly 
zeropoint drift is linear, GPC ties together all calibration fields and all survey nights 
simultaneously, producing a globally optimized solution. Calibration solutions for the 
Northern and Southern hemisphere observatories are found separately, and are tested 
for global consistency based on common fields near the celestial equator. 

Several results from the GPC are presented, including establishing candidate 
secondary standards, monitoring of near- infrared atmospheric extinction coefficients, 
and verification of global validity of the standards. The solution gives long-term 
averages of the atmospheric extinction coefficients, Aj = 0.096, Ah = 0.026, 
Ak, = 0.066 (North) and Aj = 0.092, Ah = 0.031, Ak, = 0.065 (South), with 
formal error of 0.001. The residuals show small seasonal variations, most likely due to 
changing atmospheric content of water vapor. Extension of the GPC to ~ 100 field 
stars in each of the 35 calibration fields yields a catalog of more than two thousand 
photometric standards ranging from 10*^ to 14*^^ magnitude, with photometry that is 
globally consistent to ~ 1%. 
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1. Introduction 

The Two Micron All Sky Survey (2MASS) maps the entire sky in three near-infrared (NIR) 
bands, J (1.25/im), H (1.65/im), and (2.16/um). The survey operates in both hemispheres, 
with the Northern facility at Mt. Hopkins, AZ and the Southern site at Cerro Tololo, Chile. 
Both sites acquire data with nearly identical 1.3 m cassegrain equatorial telescopes optimized 
for efficient sky coverage. Each telescope is equipped with three NICM0S3 arrays capable of 
simultaneous observations in three NIR bands. The Northern 2MASS facility began routine survey 
observations in June 1997, while the Southern facility started in March 1998. 

Nightly photometric calibration for 2MASS is derived using repeated observations of 1° x 8.5' 
calibration "fields" (see Table Q). Each hour, one of 40 calibration fields is scanned six times 
in the normal survey mode, providing six independent measurements of the standard stars in 
each field. The 2MASS calibration fields were selected initially to contain at least one "primary" 
calibration star drawn from the list of faint infrared standards published by Persson et al. (1998) 
and Casali &; Hawarden (1992). One of the fields, 92049, does not contain an a priori standard 
star, but was selected to fill a gap in the right ascension coverage of the fields. Photometry for the 
calibration stars in this field was developed using the techniques described in §3.1 of this work. 
Table |l] contains a listing of the 2MASS calibration fields and the a priori catalog magnitudes for 
the primary calibration stars in each. The fields centers listed are the coordinates of the primary 
calibration stars. The analysis described in this paper uses 2MASS pipeline source extractions 
and photometry from repeated observations of 35 of the calibration fields. 

The paper has three major parts: fQ, which introduces the global photometric calibration 
method and compares it to nightly calibration; which presents the results of major applications 
of the procedure to 2MASS data; and which restates and summarizes the main results of the 
analysis. The major applications of GPC are: (i) establishment of secondary standards (§p.l|), 
(ii) study of atmospheric extinction in near- infrared (j j3.2| ), and (iii) analysis of global consistency 
(§ |3.3D . Technical aspects of the method are given in Appendix^. 



2. Global Photometric Calibration (GPC) 

The purpose of nightly photometric calibration is to transfer instrumental magnitudes, niinst, 
onto a uniform photometric system rricai- 

rricai = minst + Ci - C2{X - 1), (1) 

where Ci is the photometric zero point, C2 is the atmospheric extinction correction, and X is 
the airmass. Figure |l] plots a typical night of calibration data, showing both the data points 



^The Kg band (pronounced 'K-short') is described in e.g., Persson et al. (1998). 
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Table 1: 2MASS calibration fields. Each field is 8.5' x 1°, centered on given coordinates of the 
fiducial standard. The equatorial fields are observed by both 2MASS facilities in the course of the 
survey to ensure uniformity between hemispheres. The listed magnitudes are the ones contained in 
the literature. 
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=*NICMOS Standards (Persson et al. 1998); 
•^UKIRT Faint Standards (CasaH & Hawarden 1992); 
■^Faint Red Standards (Persson et al. 1998); 

■^Based on 2MASS observations, red standard L547 was found variable and was subsequently replaced by 



another star, selected from secondary standards, see §3.1 



"Internal calibration field defined by the 2MASS Project; contains no fiducial standard. 



- 4 - 



(photometry of primary standards) and the nightly cahbration solution Ci (straight line). 
Each cluster of points in the nightly calibration plot is derived from one set of calibration 
observations. Each cluster consists of six data points which plot the zeropoint offset (i.e., 
meal — T^inst + C2 (X — 1)) derived for each of the six apparitions of the primary calibration star. 
Because the Survey's 2.0" camera pixels are larger than the typical ^1" seeing disk, pixelization 
flux errors dominate the scatter of the bright primary standard observations and limit the rms 
uncertainty for multiple apparitions of bright stars in 2MASS observations to ~ 0.02 — 0.03 
magnitude - underscoring the advantage of using several flux calibrators in each field to minimize 
the zeropoint uncertainty of a single calibration observation. 

Experience has shown that the nightly photometric zero point Ci is usually either a constant, 
or a linear function of time, Ci{t) = A + Bt. Equation (||) is written for a set of primary calibrators 
with known nicai (Table ||), independently for each NIR band. Solving the calibration equation (|l|) 
by a standard least squares algorithm produces a nightly photometric calibration solution, i.e. Ci 
and C2. 

The main difference between nightly and global calibration procedures is in the scope of the 
data upon which the solution is based: the former derives the calibration parameters Ci and C2 
from a single night of data, while the latter uses all survey nights. The calibration equation ([l|) is 
the same in both cases, but for GPC it can be rewritten to underscore its global nature: 

mfjki = <k + an + hn^tjn + A {Xjn " 1) , (2) 

where mf-^^ is the instrumental magnitude of the star in the A;*^ calibration field, observed 
at time moment tj during the n*'* survey night, m?^' is the true photometric magnitude of the 
star, an and 6„ are nightly photometric constants (one pair for each night), and Atjn is the time 
offset between time tj and midnight. The atmospheric extinction coefficient is now denoted by A. 
Equation (^) implicitly assumes that atmospheric extinction is constant. Seasonal variations in 



extinction coefficients (see §3.2) are thereby smoothed and represented by their average values. 
The equation is applied to each photometric band separately and produces one solution for each 
band. Note that the global procedure does not rely on a priori knowledge of rricai- the procedure 
can solve for these magnitudes as well. However, in that case the calibrated magnitudes will be 
determined to within an arbitrary constant offset. 

The global calibration procedure (|2|) provides a better photometric solution than nightly 
calibration, since it minimizes residuals globally. The benefit of global calibration is conveyed by 
Figure §: each night's photometry is calibrated through nightly parameters a and b, tying fields 
observed that night to each other, and different nights are tied together by global parameters m'^"'^ 
and A. For comparison. Figure ^ shows the actual survey sky coverage of calibration fields as a 
function of time. 

In practice, the set of equations (|2|) is written for a group of primary standards and so-called 
'tracer' stars, selected among relatively bright field stars. The same number of tracers fi is used 
in each calibration field. In principle, any number of tracers can be used to derive the global 
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Fig. 1. — Nightly calibration solutions in three bands for a typical night. The x-axis units are hours. 
Data points represent calibrated photometry of primary standards, mcai—'m-inst+C2 (X— 1), straight 
lines show photometric zero point Ci (see text). Calibration fields are indicated. Note the second 
order effect in J band (non-linear behavior of sensitivity) , which is not addressed by the procedure 
but tends to be small. 
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Fig. 2. — The essence of global calibration: photometry is bootstrapped for all fields, sources and 
nights (see text). 




Fig. 3. — 2MASS calibration coverage as the function of survey day (survey day '0' is defined to be 
March 1, 1997 UT). The R.A. coordinates of the calibration fields are shown as they are observed 
in the course of the survey (cf . Figure |2|) . Crosses and circles denote calibration fields observed by 
Northern and Southern 2MASS facilities, respectively. 
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calibration solution. However, the large data sets associated with the problem prohibit using more 
than a few tracers per calibration field. The set of equations is solved by a standard linear least 
squares algorithm (see Appendix^), producing a set of photometric nightly constants a„, 6„ for 
all survey nights, values of atmospheric extinction A, and calibrated magnitudes m'^"'^ for primary 
standards and tracer stars. Once photometry is established for the primary and tracer stars in the 
field, it can be directly applied to any of the other stars in the field. 



3. Applications of GPC 

This section lists several applications of the global photometric calibration: establishment of 
secondary standards, study of atmospheric extinction and demonstration of temporal and spatial 
consistency. The discussion in this section is based on a global calibration solutions for Northern 
data taken from May 1997 to May 1999, and Southern data from March 1998 to February 
1999. Calibration solutions for North and South are calculated separately and compared through 
common calibration fields. The global solution was obtained in two steps. First, we obtained a 
preliminary calibration solution based on the entire set of nights. This preliminary solution was 
used to create nightly calibration plots similar to the one shown in Figure ||. The nightly plots 
were examined and the few with maximal deviations > 0.05 mag from the linear behavior were 
removed from the data set. We then reran the procedure on the new 'clean' data to obtain the 
final calibration solution. Northern data consisted of a total of 14,833 calibration scans of 26 fields 
taken on 298 nights. Solving equation (||) by a linear least squares algorithm amounts to inverting 
a sparse 37758 x 675 matrix. The dimensions of the matrix, N x M, are the total number of 
observations and the number of free parameters, M = + 1) + 2n + 1, where k is the number 
of calibration fields, // is the number of tracer stars in each field, and n is the number of survey 
nights. In the Southern data, there are 227 observing nights for 26 fields (11945 total calibration 
scans) and the solution is found by inverting a smaller, 33407 x 533 matrix. Both hemispheric 
solutions are obtained based on the photometry of fiducial tracers and two tracer stars in each 
calibration field (i.e., /i = 2). In solving equation (^, we look for an unconstrained solution (i.e. 
calibrated magnitudes of primary standards are free parameters; see j|A[), to ensure no bias enters 
from the known magnitudes of primary standards. Because the solution is unconstrained, the 
GPC calibrated magnitude system will have an arbitrary constant offset in each band. The offset 
for each band is found from averaging the difference between GPC solution and a priori catalog 
magnitudes of the fiducial standards in Table ||. The GPC system magnitudes for the primary 
standards in each field are given in Table g 



3.1. Secondary Standards 

To establish candidate secondary standards in the calibration fields, we apply the GPC 
solution to additional field stars. In order to exclude variable stars, secondary standard candidates 
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must have a low rms uncertainty, 



1 ^ 



\ i=l 

where N is the number of individual observations (scans) of a source. To ensure the statistical 
significance, we keep only the candidates with N > 100. Figure |^ shows the pooled rms plots for 
all solutions. Each point in the plot represents a single star. Dashed lines, which indicate the 
adopted thresholds for rms error a, were obtained by fitting a fourth-order polynomial^ fifn) to 
the data and adjusting its intercept point so that the number of points below the curve is 66% 
(Icr) of the total. Toward fainter magnitudes f{m) rises sharply, so we modify our cj-threshold to 
pick sources with 

a < min {/(m), 0.05}. 

The above selection criteria produce a sample of low-variance stars, which are good candidates for 
secondary standards. The high-variance sources that fail our cr-criterion are often variable stars, 
which can be reliably identified by examining their light curves. 

After applying the GPC solution to each hemisphere, we obtain two lists of candidate 
secondaries, one for each hemisphere. The magnitudes of stars in equatorial fields observed from 
both hemispheres are merged using variance weighted average: 

m = (4) 

where rriN and ms are average magnitudes for the star in the Northern and Southern solutions, 
and (Tat, 175 are respective root variances. The list of 2177 candidate secondary standards for all 
calibration fields is given in Table ^ in Appendix ^. This list is used by the 2MASS Project to 
select a subset of unconfused point sources for actual calibration of 2MASS data. 



3.2. Atmospheric Extinction Coefficients 

The GPC solution provides an estimate of the atmospheric extinction in all three bands. The 
atmospheric extinction in the near- infrared is caused primarily by molecular bands of H2O (Hall 
& Genet 1982, Mountain et al. 1985) and, to a lesser extent, by CO2 (Manduca & Bell 1979). In 
equation @) extinction enters through the coefficient A, which represents the average value of the 
atmospheric extinction over the period sampled (one number for each band). To decouple the 



*A third-order polynomial was used in Ks band 
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Fig. 4. — Pooled rms error plots for both hemispheric solutions. Toy row: Northern solution; Bottom row: Southern solution. 
Each dot represents one star. Dashed lines show the adopted thresholds (see text). 
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atmospheric extinction parameters from photometric zero points (sensitivities), the observations 
during each night span a range of air masses (X). A more detailed model for atmospheric 
extinction, for example involving a periodic function to mimic seasonal variations, could also 
be considered in equation (^). However, any extinction model can be devised and applied a 
posteriori, once we compute the residuals between the regressed data and the observations. 



Am = rrhf-ai — rriinst- One such model is presented in §3.3 



The atmospheric extinction coefficients for North and South produced by GPC are given 
in Table |2|, which also lists CTIO values from Frogel (1995). The uncertainties in extinction 
coefficients are random errors, obtained by inversion of the Hessian matrix. The values listed in 
the table are consistent at the 2 — a level. Figure ^ shows the residuals rriinst — meal — a — 6At as a 
function of airmass. The slopes of the linear regression lines are the mean extinction coefficients. 

Given the time baseline of 2MASS and the homogeneity and accuracy of 2MASS photometric 
data, we have a great opportunity for studying atmospheric extinction in the near-infrared over 
the period of more than a year. Figure ^ shows the seasonal variations in atmospheric extinction. 
The monthly averages were derived by, first, calculating the residuals without the extinction term 
for each month separately, and then, by fitting a straight line to the distribution of residuals 
as a function of the airmass. While some scatter is present in Figure ^, the overall behavior 
of monthly extinction averages is in agreement with expectations. The water vapor content in 
the atmosphere peaks during Northern spring and summer, which leads to higher atmospheric 
extinction in the North during that season. The amplitude of the seasonal changes is about 0.02 
magnitudes/airmass. 



3.3. Global Consistency 

The global consistency of the solution can be assessed from the analysis of residuals as a 
function of spatial coordinate (spatial uniformity) or time (temporal uniformity). In this section, 



Table 2: Atmospheric extinction coefficients and mean squared norms of the residual vectors from 
GPC. The extinction coefficients are in magnitudes per unit airmass. The CTIO values from Frogel 
(1998) are listed for comparison. 



Band 


Hemisphere 


2MASS Ax 


Frogel ^A, CTIO 


J 


North 


0.096 ± 0.001 




J 


South 


0.092 ± 0.001 


0.100 ±0.024 


H 


North 


0.026 ± 0.001 




H 


South 


0.031 ± 0.001 


0.055 ±0.021 


Ks 


North 


0.066 ± 0.001 




Ks 


South 


0.065 ± 0.001 


0.085 ±0.018 



0.5 1 0.5 1 

(X-1) (X-1) 

Fig. 5. — Density plots of the residuals (without the extinction term) as a function of the airmass. 
Left: Northern solution; Right: Southern solution. Vertical axis is in units of magnitudes. Contour 
levels are logarithmic, corresponding to point density of 0.3, 1, 3, 10 and 30 mag~^ airmass"^. 
The slopes of the straight lines are the derived mean atmospheric extinction coefficients listed in 
Table |. 

we present the results of two such analyses. First, we examine the uniformity of the standards 
between the hemispheres, from looking at residuals Am = rrif^orth — mSouth for common calibration 
fields. Then, we test spatial and temporal uniformity by comparing the calibrated magnitudes of 
primary standards (derived from GPC) with literature photometry. 



3.3.1. North vs. South Uniformity 

The linearity of the GPC solution is demonstrated in Figure |^, which shows the difference 
between catalog and GPC magnitudes for primary standards as a function of the catalog 
magnitude. The figure shows that the absolute value of the difference for most primaries is 
less than 0.02 magnitude, i.e. consistent with zero. The absence of any trends in these data 
demonstrates the linearity of the GPC solution. Note that the results of such comparison are 
limited by the accuracy of input catalog since the errors in the fiducial magnitudes of primary 
standards are propagated unchanged. In fact, the good agreement demonstrates the high accuracy 
of the input catalog. 

The comparison of the photometry between the hemispheres is carried out in equatorial 
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Fig. 6. — Seasonal variations in the atmospheric extinction (in mag/airmass) in the North [left 
panel) and in the South {right panel). Monthly average coefficients in J, H, and Kg are plotted 
with triangles, squares, and circles, respectively. The errors are indicated. Horizontal dashed lines 
show global averages in three bands derived from the respective GPC solutions. No August data 
for the Northern solution were available. 

fields (see Table |^). Data from each hemisphere were calibrated separately by unconstrained 
GPC. The solutions (nightly photometric constants a and b, and extinction parameter A) were 
applied to the 100 brightest stars in each calibration field. Figure ^ shows the difference between 
calibrated magnitudes of common stars as the function of coordinate. Except for a few outliers, 
the largest excursions for individual secondary standards (crosses) are bounded by Am < 0.04, 
while the largest excursions for the primary standards (squares) are less than 0.02 mag. The 
mean differences for the photometry of common standards (squares in Figure P) are 0.002 ± 0.007 
(J), -0.002 ± 0.007 (H), and -0.001 ± 0.008 {K^). The plot of North-South photometric offset 
vs. declination displays a slight trend of unknown origin. The slopes of the corresponding 
regression lines in J, H, and Ks are (—4.9 it 0.9) x 10~^ mag/deg, (—7.5 ± 0.9) x 10"^ mag/deg, 
and (—3.9 ± 0.6) x 10~^ mag/deg, respectively. The same regression analysis using primary 
standards only gives slopes consistent with zero at the 2a level: the corresponding slopes are 
(-4.9 ±2.8) x 10"^ in J, (-5.2 ±2.8) x 10"^ in H, and (-0.4 ±3.5) x lO""^ in K^. Figure | shows 
the same residuals as in Figure ^, except as a function of magnitude in the North. The slopes of 
the regression lines in all three bands are consistent with zero. 
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Fig. 7. — The difference between catalog and GPC magnitudes as a function of catalog magnitude 
for primary standards. Symbols represent J, H, and Kg magnitudes (crosses, triangles and squares, 
respectively). Most data points differ from zero by less than 0.02 mag. Note the absence of any 
trends over the entire magnitude range. 



3.3.2. Temporal Stability 

To analyze the temporal stability, we consider the difference between calibrated and fiducial 
magnitudes of primary standards as a function of time. Specifically, we calculate the difference 
'calibrated-fiducial' for all primary standards observed on a given night and plot the average 



difference. Figure 10 shows the corresponding differences in all three bands for Northern and 



Southern hemispheres. The figure, which essentially indicates the long-term stability of the global 
calibration, shows a very small temporal variation in the photometry of the primary standards. 
The amplitude of the variations is less than 1%. 



4. Summary 

• We have presented a global photometric calibration procedure and applied it to more than 
a year of 2MASS calibration data. Exploiting the fact that the observed nightly zeropoint 
drift is linear, one can simultaneously calibrate the photometry of fiducial standards and any 
number of tracer stars. The solution is found by using a standard least squares algorithm 
and is globally optimized, in the sense that it is best solution for all survey fields and nights. 
The solution does not presume a priori knowledge of the magnitudes for any of the stars, 
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Fig. 8. — 'North-South' difference between calibrated magnitudes of stars in common fields as a 
function of right ascension {left panel) and declination (right panel). Crosses represent field stars, 
squares represent primary standards. Clusters of points are calibration fields. Due to elongated 
shape of calibration fields in latitude, the clusters are more dispersed in the right panel. Dotted 
horizontal lines are drawn at = ±0.04, dashed line is at zero. 



but instead produces the optimal set of relative magnitudes. A constant offset is determined 
by averaging the differences between the calibrated and catalog magnitudes for the primary 
standards. 

• The photometry of the primary and two tracer stars per field is transferred to the 100 
brightest non-variable stars in each field. With a u-threshold, we select a low-dispersion 
sample of 2177 candidate secondary standards to accompany the set of primary 2MASS 
standards. The selected secondaries fall in the magnitude range 9.0 < J < 14.5, 

9.0 <H < 13.9, 9.0 <Ks< 13.7. 

• Using the global solution, we derive the atmospheric extinction. The global averages are 
Aj/Ah/Ak, = 0.096/0.026/0.066 (North), and Aj/Ah/Ak, = 0.092/0.031/0.065 (South). 
The averages are consistent with previous work at the 2a level. 

• Seasonal variations in the atmospheric extinction are analyzed. The amplitude of the 
variations is approximately 0.02 mag/airmass and is consistent with the behavior expected 
from variations of the water vapor content in the atmosphere. 




Fig. 9. — Difference (in magnitudes) between Northern and Southern solutions as a function of 
Northern magnitude. The notations are the same as in Figure |^. 
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Fig. 10. — Difference between calibrated and fiducial magnitude for primary standards as a function 
of time for Northern (left panel) and Southern (right panel) solutions. Each survey night is 
represented by a single point which displays the average difference for all standards observed on 
that night. 
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A. Least Squares 



The set of equations @ represents a hnear least squares (LLS) problem. The general LLS 
problem is well-documented in the literature (e.g., Lawson & Hanson 1973, Greenbaum 1986) 
and has many professional codes available for its solution. Written in the matrix form, the 
mathematical model for the GPC is 

^>x = b, (Al) 

where ^ is the design matrix, x is the parameter vector and b is the vector of observed 
(instrumental) magnitudes. The parameter vector x consists of the magnitudes of the fiducial 
standards in each field, followed by the magnitudes of the tracer stars for each field, followed by 
the photometric nightly constants a, b for each night and the atmospheric extinction coefficient A: 



(mf^...,m|*^m?i. 



> bn^ -A 



(A2) 



where k denotes the total number of calibration fields and n denotes the total number of 
observation nights. The number of tracers in the field is fi (the same number of tracers in each 
field is used). The design matrix $ has the following structure: 
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The tracer stars are grouped by the fields, i.e. after the fiducial standards, there are n tracers 
from the first field, then there are n tracers from the second field, and so on. There is a single 
unity in the first k{p, + 1) columns of each row of the design matrix, representing the observation 
of either the fiducial standard or a tracer star. The subscript j denotes the time offset from the 
midnight on the night re. The dimensions of the matrix are N x M, where N is the total number 
of observations of all fiducial standards and tracers, and M = A;(/x + 1) + 2n + 1 is the number of 
free parameters. 



This preprint was prepared with the AAS lATp^X macros v4.0. 
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A.l. Constrained Scheme 

The least squares solution of the equations is globally optimal (i.e., it minimizes rms), 
but, generally speaking, has an arbitrary offset from the true magnitude: 

A,o\ _ / A,0 



To produce the optimal solution without the arbitrary additive constant Kn. one has to set the 
magnitudes of the fiducial standards as fixed points of the GPC. This results into GPC with 
constraints, where the constraint equations are: 

r^L = ^fk'' (A4) 

(the superscript std denotes the fiducial standards in the fields). In the matrix form, the calibration 
problem with constraints is written as 

<l>x = b, 
Cx = d. 



For this particular problem, the constraint equations (A4) result in the identity constraint matrix 
C, and the LLS problem with constraints is straightforward to solve using, e.g., orthogonal basis 
for the null space of the constraint matrix (Lawson & Hanson 1973, Ch. 20). The GPC with fixed 
points derives the globally optimal solution which is the closest to the listed magnitudes of the 
fiducial standards. 



B. Secondary Standards 



Below, we present a sample table of 2MASS secondary standards, which includes the primary 
and some of the secondaries for the field 90021. The primary standard is the first star listed in 
each field in Table ^. Columns in the table are aj2ooo, (572000 coordinates of the stars, the globally 
calibrated JHKg magnitudes with the corresponding uncertainties, and the number of observations 
(i.e. individual scans) in each band. The positions of the secondaries are accurate with respect to the 
ICRS to < 0.5" rms. Finding charts for the secondaries can be obtained by using 2MASS Visualizer 



tool at [http : / / ir sat est . ipac . caltech. edu: 8001/ applications/2MASS/ReleaseVis/|\ 
The full table of candidate secondary standards in all calibration fields is 
available from FTP archives at ftp : //nova, astro .umass . edu/pub/nikolaev/ , or at 



ftp : //anon-f tp . ipac . caltech. edu/pub/2mass/globalcal/ 



^Finding charts are available only for calibration fields which are covered by the Second Incremental Data Release 
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Table 3: 2MASS Secondary Standards (sample). Columns list J2000.0 coordinates of the stars. 
JHKg magnitudes and their rms errors and the number of observations in each band. 
a (J2000) 5 (J2000) J H ax, Nj Ng Nk, 

Field: 90021 
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